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Abstract 



Using a' 3 terms of type IIB supergravity action we study higher order corrections to the non- 
relativistic non-extremal D3-brane. Utilizing the corrected solution we evaluate corrections to 
temperature, entropy and shear viscosity. We also compute the rj/s ratio which although within 
the range of validity of the supergravity approximation and in the lowest order of the correction 
the universal bound is respected, there is a possibility for a violation of the bound when higher 
terms in the expansion are taken into account. 



1 Introduction 



AdS / CFT correspondence [1] has been found to be a powerful tool for studying strongly coupled 
field theories in terms of weakly coupled gravities. For instance it has provided us a framework to 
study thermodynamic and hydrodynamic properties of certain strongly coupled gauge theories 
which have gravity dual (see for example [2]). In particular in this context the AdS/CFT 
correspondence has been used to evaluate the ratio of shear viscosity to the entropy density 
(see for example [3-12]). 

For example, for M = 4 supersymmetric Yang-Mills theory in 4D at finite temperature 
whose gravity dual is given in terms of the near horizon of the non-extremal D3-brane it was 
shown in [13] that at large t'Hooft coupling and for large number of colors (large N) the gravity 
calculations give 7 = ^ which is compatible with the RHIC data [14]. In fact motivated by 
gravity calculations it was proposed [3] that the rj/s is bounded from below 

S 47T 

Further investigations have confirmed this, rather universal, behavior for broader examples [15]. 
We note, however, that since in the gravity side we have higher derivative corrections to the 
tree level action which has been used to reach to the above conclusion, one may wonder if 
these corrections would spoil the bound. Indeed this point has been studied in several papers 
including [16-19]) where it was shown that there is a possibility for the bound to be violated 
when higher order corrections are taken in to account. 

In this article we would like to extend the above considerations for those non-relativistic 
field theories whose gravity duals are given in [20,21]. The corresponding gravity solutions may 
be embedded in ten dimensional type II super gravities. Indeed starting from brane solutions 
in type II super gravities and using a Null Melvin Twist [22, 23] one can find new solutions 
(non-relativistic branes) which provide gravity duals for non-relativistic field theories [24-26]. 
Following the relativistic case one would expect that heating up the non-relativistic field theory 
corresponds to adding a black hole solution to the bulk [27-30]. 

Having had the gravity description of the non-relativistic field theory one may proceed to 
study thermodynamic and hydrodynamic properties of the theory. Indeed at leading order this 
has been done in [31,32] where it was shown that the ratio of the shear viscosity to entropy 
density is the same as that in the relativistic case; f]/s = I /Air. It is then natural to pose the 
question whether there is a lower bound for this ratio in the non-relativistic case as well. To 
explore this question we will study the effects of higher order corrections to this ratio in the 
context of AdS/NRCFT. 

To be specific we will consider non-relativistic non-extremal D3-brane in type IIB string 
theory which can be used to study thermodynamic and hydrodynamic properties of the dual 
three dimensional non-relativistic field theory. In leading order the gravity solution can be 
obtained from D3-brane solution by making use of the Null Melvin Twist. The resultant 
background is still a solution to the leading order action of type IIB supergravity. Now the 
task is to find the corrections to the solution when a' 3 terms are added to the action. This can 



1 



be used to study the effects of higher order terms to the 77/s ratio 

We note, however, that in general it is difficult to find a closed form for the corrected 
solution. Nevertheless one can solve the equations of motion in the presence of higher derivative 
terms perturbatively. In fact there are two parameters that control the perturbation. The first 
one, denoted by 7, controls the higher order terms in the action. The other, denoted by A, 
parametrizes the deviation of the geometry from that in the relativistic case 2 . Taking these two 
as the expansion parameters one finds the following schematic expansion for a typical quantity 
Q 



Q=Qo 



1+ K + cl (\ 2 T)+c( (A^T + .-. 7 + - 



(1.2) 



120 - 183.31tt 2 (A 2 T) 2 + 0((\ 2 TY 
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(1.3) 



where are some numerical factors. In particular for t]/s ratio at first subleading order we 
find 

rj 1 

s 4ir 

which has a potential to violate the bound. We will back to this point later. 

The paper is organized as follows. In the next section we will review the tree level non- 
relativistic non-extremal D3-brane to fix our notation. Then we will find higher order cor- 
rections to the solution when the a' 3 terms are added. In section three we will study the 
thermodynamic and hydrodynamic properties of the dual non-relativistic theory in the pres- 
ence of higher derivative terms. In particular we will find the corrections to the temperature, 
entropy and shear viscosity. The last section is devoted to conclusions where we also give 
comments on the rj/s ratio. 



2 Higher derivative corrections 

In this section we will study higher derivative corrections to the non-relativistic non-extremal 
D3-brane. To proceed we will first fix our notation by reviewing the tree level action of type 
IIB supergravity and its non-extremal D3-brane solution. When only dilaton, metric, B-field 
and RR-five form are non-zero the tree level action of type IIB supergravity in the string frame 
is given by 

I = r-V / d 10 xy^gC 

= 4/ A ^{ e "( fl+4a, '* s ''*-^)-^ F5 »}- (2 - l) 

The near horizon limit of the non-extremal D3-brane solution may be written as follows 

ds 2 = r 2 (-fl(r)dt 2 + dy 2 + dx\ + dxj) + - + (dx + A) 2 + ds 2 cp2 , e" 2 * = 1 , 

Jo \ r ) r 

1 Higher order corrections to the non-relativistic background have recently been studied in [33] where the 
higher derivative terms were given by Gauss-Bonnet action. It was then shown that the higher order terms may 
correct the power exponent and as a result the rj/s bound may be violated. 

2 To be precise the relevant expansion parameter is A 2 T where T is the temperature. 
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F {5) = dC {A) = 2(1 + *)d X A J A J, dA = 2J, V ol(C P 2 ) = \j A J . (2.2) 

Here /q (r) = 1 — r 4 /r 4 and 

A = - sin 2 yu(c?^ + cos 6d<p) , ds 2 cp2 = dfi 2 + — sin 2 /i(a 2 + a 2 + cos 2 /iO" 2 ) , (2.3) 

where Oi are the SU(2) left invariant one-form 

Ci = cos ipdO + sin 6 sin ^<i</> , cr 2 = — sin ipdO + sin 6 cos ?/?<i(/> , as = dtp + cos 9d(fi . (2.4) 

Using the Null Melvin Twist one can map the solution (2.2) to the non-relativistic non- 
extremal D3-brane which is still a solution of the tree level type IIB supergravity action. The 
obtained solution is [25] 

dg 2 = r^fpr) / + r2x2)df2 + 1 - ^ 2 f° (r V - 2r 2 A 2 ^ 



k(r) V V ' f§{r) 

dr 2 1 

+ ^(dxl + dx^ + j^^—idx + Af + ds 2 ^, 

B {2) = ^(d X + A)A(f (r)dt + dy), e 2 * = ^ , (2.5) 

where k(r) = 1 + r 2 A 2 (l — /o (r)) an d the RR-five form remains unchanged. 

Now the aim is to study the effects of higher derivative terms to the non-relativistic non- 
extremal D3-brane solution. Here we restrict ourselves to corrections from gravity side. There 
are other efforts have been done to include other correction terms, for example see [17]. The 
higher derivative terms which we are interested in are the a' 3 corrections to type IIB super- 
gravity action. In string frame they are given by 

h = 7TT f d w x^ 1 C 1 = ttV / ^^7^, (2.6) 



where 7 = |C(3)(«') 3 and W can be written in terms of the Weyl tensors 

VV — Lj ^pmnq^h ^ rsk ~T~ ^ ^pqm.n^h ^ rsk ■ 

In principle one needs to solve the equations of motion coming from the action Iq + I\. Of 
course, in general, it is not an easy task to solve the resultant equations. Nevertheless we may 
start from an ansatz and try to solve the equations of motion for the parameters of the ansatz. 
We note, however, that in general the obtained equations cannot be solved exactly, though they 
may be solved perturbatively. 

For making our ansatz, we start from the same general ansatz as [34] for relativistic non- 
extremal D3-branes. Then we do all the steps for Null Melving Twisting of _D p -branes noted 
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in [25] for this general background. Finally we find the following ansatz for the non-relativistic 
non-extremal D3-branes 



ds 2 = -(1 + r 2 X 2 s 2 (r))dt 2 + "\ JJ — dy 2 - 2r 2 X 2 s 2 (r)dtdy 

k(r) V f (r) 



f ^-P-dr 2 + S -^l(d X + A) 2 + s 2 {r)ds 2 CP , 
K{r) 

e 2v(r) 



2<i> 



F {5) = dC {4) = 2(1 + *)d X A J A J , e 

= ^^(dx + A)A(/(r)<ft + dy), (2.8) 

for arbitrary functions f(r), h(r), s(r) and u(r). Moreover 

jfe(r) = l + r 2 AV(r)(l- f 2 (r)). (2.9) 

Note that in the limit of A — > the above ansatz reduces to that in the relativistic case. We 
will also assume that f(r) and h(r) have zero at r = ±ro and the extremal limit is obtained at 

Using the above ansatz the equations of motion of the total action Jo + I\ read 

dC d dC d 2 dC _ dd d d£ x d 2 dC x 

dXi dr dX[ dr 2 dX'l ~ 1[ 8Xi dr dX[ dr 2 dX'/' ' l ' } 

where prime denotes the derivative with respect to r and Xj = {f(r),h(r),s(r),v(r)}. To 
proceed we choose the following perturbed ansatz over the tree level solution 

f(r) = (l-^)l(l + 7 F(r)), /i(r) = (l-^)-4 (1 + 7^(0) , 

s (r) = 1 + 7%), u(r)=7^(r). (2.11) 

It is straightforward, though messy, to plug the ansatz to the equations of motion to find the 
following differential equations for the perturbed functions {F(r), H(r), S(r),V(r)}. Setting 
p = — and A = Ar one gets 



'■0 



(P 4 - 1) (PS" - \ P V" - + 4 (p 4 - i) (S> - \v>) - 4p 3 {s-\v+ \h) 

F ^ (2.12) 



p 13 (p 2 + A 2 ) 5 



p (p 4 - 1) (F" + 45" - 2V") + 10(p 4 - i) (S' - 2V) - 4 (p 4 - i) H' 
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+ P H6 P F' + 52S-SV-20H) = -^M^ , (2 .13) 

(P 4 " 1) r + ^ - 1) (55' - 2V0 - V (H + 5 S - \V) = , (2.14) 

p (p 4 - 1) (F" + 5S" - 2V") + 5(p" - i) (55' - 2V') - 4(p 4 - i) 

A 

+ p 3 (6pF' - 20H + 20S) = A) . (2.15) 

p 13 (p 2 + A 2 ) 8 

Here Fi are functions of p whose explicit forms are given in appendix A. Now the task is to 
perturbatively solve the above equations with given boundary conditions. First of all we would 
like to have a solution whose asymptotic is the same as that for the non-relativistic non-extremal 
D3-brane. Moreover we require that the solution to be finite at the horizon. 

For our purpose we just need to solve the equations up to (9(A 2 ). 3 The results are 



640 - 1375p 4 - 1465p 8 / 11978, 4 2837, , ,„ 1, 6304 
32^ " + (— <" " 5989> ln(1 + ?» " IT 



1564664 11978 , 94 1 2197 1 109523 1 2332 1 \ A 



1575 15 P+ 45p 2 + 450 p 6 350 p 10 + 21 p 14 ) p 4 - 1' 



2480p 8 + 1895p 4 - 10000 /,5989 „ 9036 4 2837, , , 1 , 
= 64/* + ( ( — " " + — )ln(1 + ^ 

6304 , , , 1657358 18777 1 45493 1 607097 1 24803 1 
+ — — ln(2) 



15 v ' 1575 28 p 14 75 p 10 6300 p 6 75 p 2 

21119 11978 5989 \X 
15 H 25 ' 5 H ) p 4 - 1 V ; 



15 + 90p 4 /,5989 , 3257, , , l x 5989 , 3782 1 1532 1 



64p 8 V 25 25 P 25 75 P 2 375 P 

184273 1 53 1 

10500 7° ~ ^P 11 



A 2 . (2.18] 



= 30 + 45/ + 90/ 1 / + + J_ _ 735Q J_ _ 24M J_ 

16p iz: 175 V p^ p z p D 



5 Note that using the expression of the temperature, (3.2), at leading order one has A = Ar = A 2 T. 



5 



- 1617o4tt + 3375^7 M 2 • (2-19) 

These are all information we need to evaluate the corrections to the temperature, entropy and 
shear viscosity which we will do in the next section. It is worth noting that the above solution, 
after going to Einstein frame and using a field redefinition, reduces to the known results in [34] 
for A = . 



3 Entropy and Viscosity 

In this section we will study thermodynamic and hydrodynamic properties of the finite temper- 
ature non-relativistic three dimensional field theory whose gravity dual is given in terms of the 
non-relativistic non-extremal D3-brane. In fact, by making use of the results of the previous 
section, we will be able to study the effects of the higher derivative corrections to the entropy 
and the shear viscosity. To proceed let us start from the thermodynamical quantities. 



Temperature 

Since in the presence of the higher derivative terms the solution gets correction, one expects that 
the temperature will also be corrected. To find the temperature one may follow the standard 
route. Namely we will have to identify the Euclidean time in the corrected geometry in such 
a way that the corrected metric will be free of conical singularity. Consider as the null 
generator of the horizon, the surface gravity is defined by 



K 



(3.1) 



and the temperature is related to surface gravity by T = As indicated in [27] the Killing 
generator of the horizon is given by £ = -^J^ which differs from the relativistic case by a 
prefactor. Knowing these, to first order of 7 and A 2 one finds 



T 



1225 
(19.14 



7) 



,3152 



v 15 
14.35A 2 )7 



ln(2) 



206806 
1575 



) 7 A 2 



ttA 



vrA 



(3.2) 



Note that in the limit of 7 — > this reduces to known leading order expression [27]. On the 
other hand for A — > and AT — ► T we recover the relativistic result as well [5,34]. 



Entropy 

Entropy is one of the interesting quantities usually people study in this context. To compute 
the entropy one may use the Euclidean method by which we need to evaluate the value of the 
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action on our solution. We note, however, that in general the action diverges and needs to be 
regularized by subtracting the value of the action evaluated on a reference background. To do 
this we will have to be careful about how to define the cut off. Although this method is found 
useful for geometries which are asymptotically AdS, in general for geometries with arbitrary 
asymptotic, this procedure is not easy to be applied. Therefore it will be useful if one can use 
another method. 

Indeed there is an alternative way to compute the entropy via Wald's formula for the 
entropy [35-37]. When the Lagrangian does not contain the covariant derivative of the curvature 
the Wald formula for the entropy reduces to the following simple form [38] 

Sbh= / dx H ^-±—g^gi p , (3.3) 

where C is the total Lagrangian and the integration is over the horizon, g is the orthogonal 
metric to the horizon which can be obtained in terms of the normal vectors to the horizon as 
follows 

g^ = {N t ),{N t ) u + {N r ),{N r ) u . (3.4) 
In our case setting £ = ^-(£ + 7^i) an d taking into account that 

N t = y/gTt(l,— ,0,0,0,0), ^. = (0,0,0,0,^,0), (3.5) 



9tt 

the Wald entropy (3.3) reads (for more technical details see e.g. [39]) 



2ni 



% n f h r~n [ ®£ ± ± , | dC , | I 

-2- / dx vV < T7p 9tt 9rr + ~pTB 9yy 9rr + JTp 9ty 9rr f ■ (3-6) 

10 J H L OrLtrtr Oity r y r Oitt r y r J 



Using the corrected geometry presented in the previous section it is easy to compute each term 
in the above expression for the entropy. Altogether up to 0{\ 2 ) one arrives at 



2n 2 RV 2 V CP 2V l ,/„ 4635 , ,12608, 1356424, T2 



2tt 2 RV 2 V CP 2V 1 3 
2^2 



'10 



rl (4 + 289.68 7) - 835.82 7 A 2 . (3.7) 



Here R is the radius of compact direction y, V± is the volume of the x coordinate and V 2 is the 
volume of xi and x 2 directions. Using the expression for the temperature (3.2) we can eliminate 
r to express the entropy in terms of the temperature 



OBH - 7T/2, l AJ ) 



1 + 15 -252tt 2 (X z TY 7 



(31 
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Shear viscosity 

The shear viscosity of a translation invariant theory at the hydro dynamic limit can be evaluated 
by making use of the Kubo formula [40] 

^ = -limIcj G !« ( W) fc = 0), (3.9) 

UJ — >\J (jj 

where is the retarded Green's function of the component T XlX2 of the energy mo- 

mentum tensor. To compute the retarded Green's function we utilize the AdS/CFT dictio- 
nary [40,41], by which the retarded Green's function can be related to the on-shell bulk action 

G* lW >,fc) = lim 2F(u, k, r) . (3.10) 



r— >oo 



The shear viscosity in the non-relativistic field theory has been evaluated at leading order 
in [32] where it was shown that the ratio of the shear viscosity to the entropy density saturates 
the universal bound 

5 = i (3.11) 
s Air 

Following [5] we would like to study the effects of higher derivative corrections to the value of 
the shear viscosity. In order to evaluate the shear viscosity we consider a perturbation to the 
metric as h XlX2 = r 2 ip(t, r, x\, x 2 ) which can be Fourier transformed along t and x^s directions 
leading to 

V(t,r,y, xi ,x 2 ) = J ^0 e - Mk ^ jk (r) . (3.12) 

Since in the Kubo formula, the spatial momentum is zero, we restrict our calculations to k — 0. 
On the other hand since the energy momentum tensor has zero particle number the perturbation 
must be y-independent. Actually it is easy to see that the relevant perturbation of the metric 
decouples from other perturbations and in fact it obeys the equation of motion of a minimally 
coupled massless scalar [16] which in our case we have 

M + dpi + 2D<p u - ^(Ftpl + 2B<pT u + C(p u ) + ^r 2 WL + + A<p u ) = . (3.13) 

The coefficients A, ...,F are given in the appendix B. We note, however, that in general it 
is difficult to find an explicit solution for the above equation. Nevertheless we may solve it 
perturbatively. Changing the variables as before, p = ^ and A = Ar , and expanding the above 

equation up to 0(A 2 ), one finds 

P(P 4 - l) 2 ^ + (5p 8 - 6p 4 + 1)^ - 4p 3 fi 2 (AV - p 2 - A 2 K 

+ Vv _ i) + ° 2 ^r + c ^ + c ^ + = » ( 3 - 14 ) 
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where the coefficients C, are given in appendix C. Next we set 



<p u (p) = (1 - -,) A G^p) 



(3.15) 



where G u {p) is a regular function at horizon and A is a constant to be determined. Setting 
uj = 2r tt>, for a solution normalized to one at p — > oo up to the order of O(vo 2 ) we have 



MP) = (l-^{l+7^(l201n(l + i)-^ 



U 120 ■ (3.16) 



p 6 p 10 ; 



1 . 7 .3385 2335 768. ~ 2tia . . 
A(») ( ln(l + -) - ^ + __ + _) + 7 A 2 lU(p) 



16 v p 4 



Here VK(p) is an unknown function satisfying the following differential equation 



dp 2 



dp 



15 



+ ~ ^96p 24 - 660p 20 + 1344p 16 + 248p 12 - 2368p 8 + 1340p 4 ^ ln(l + 



32 1 p 22 - 29p zu + 



7278127 18 9859 16 827 14 184637 12 431 10 



12600 



■P 



15 



120 



360 



120' 



12381 14 o 201 fi 96023 . 9205 12608 4 1Sl fA . 

p 8 + p 6 + p 4 p 18 n 4p) = . 

200 1 40 P 56 P 8 15 P V HJ 



(3.17) 



Using Maple one can solve the above equation which in turn can be used to find the solution 
<p ro - On the other hand the constant A can be fixed by the regularity condition on G u (p). 
Following [19] we note that in order to have the regularity up to order 7, we need to check the 
regularity condition up to order of 0(tt> 2 ). Doing so one finds 



A(») 



2ir Q vo 
AttT 



2A 



15 



3152, , n . „ 2 206806 
■ln(2)r 2 A 2 + 



1575 



64 J 



(3.18) 



Using the equation of motion (3.13) one finds the quadratic part of the on-shell action of the 
bulk theory as follows (for more detail see [5]) 



Vcp 2 Vi 
2ac 2 



F' \ 1 

B-A- — + 2p E\ ifiip^ + -(C- A')<p u (p- U + E Pl (p' u (p'_ 



(3.19) 



Putting everything together we get 



G 



X 1 X 2 ,X 1 X 2 



vo,0) = lim 2jF(to, 0,p)= lim — 



.2 

'10 



(2p 2 A 2 -l+ivo)+ p 4 
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+ 



+ 



3232 



tvo — 



2825 



51712 

Finally the shear viscosity becomes 



+ 



394 ,- 2l , . 

)A n 2 

303 ; V ; 

7325 



220583, 380537 



169680 v 661749 



+ itx>)\ 2 



2V CP 2V 1 w 4 
ri = hm 7, r n I 1 



25856 
12315 



-1 



61575 
206848 

,3232 



-ivo) 



5989 
8080 



P 4 A 2 )7 



ln(2) 



441166 



) 7 A 5 



(3.20) 



(3.21; 



64 ' v 5 y ' 525 

Using the expression for the temperature (3.2) we can eliminate tq to write the shear viscosity 
in terms of the temperature as follows 



V 



-(AT? 



135 -435.31vr 2 (A 2 T) 2 7 



(3.22) 



4 Conclusions 



In this paper we have studied the non-relativistic non-extremal D3-brane in the presence of 
higher order corrections given by a' 3 terms in type IIB supergravity. By making use of the 
corrected solution we have computed the effects of the higher order corrections to the temper- 
ature, entropy and shear viscosity. Having computed the corrected entropy and shear viscosity 
it is natural to obtain the ratio of the shear viscosity to entropy density. To find it, we first 
need the entropy density which turns out to be 

A^VcpiV! 



(AT? 



1 + 15 - 2527^(A^Tr 7 



(4.1) 



Using the expression for shear viscosity one finds 



1 
s 



1 

47T 

1 
47T 



1 + 1207+ 161n(2) 



972 



tt 2 7(A 2 T) ; 



120 - 183.31tt 2 (A 2 T) 2 7 



(4.2) 



which presents the correction to the tree level bound 1/Att. Note that in the limit of A 2 T = 
we recover the result of the relativistic field theory [19]. When 7 = we get t]/s = 1/Att which 
is independent of A showing that at tree level as far as the ratio is concerned both relativistic 
and non-relativistic cases lead to a universal result [32]. 

It is now important to see whether this correction can violate the lower bound of the rj/s 
ratio. This depends on the sign of the 7 term, 120 — 183.317T 2 (A 2 T) 2 . In fact for A 2 T > 0.257 
the coefficient of 7 is negative leading to a violation of the bound. On the other hand for 
A 2 T < 0.257 the higher order corrections respect the bound. Therefore it is crucial to estimate 
A 2 T as precise as we can. To find a bound on the value of the corrections we note that the 
supergravity approximation is applicable when the curvature is small 4 which in our notation 



4 In our notation we have set the radius of the space time to one. 
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it means \1Z\ = 20?"qA 2 <C 1, in other words, A 2 T <C 0.07. Therefore within the supergravity 
approximation and to the order of 0(X A T 2 ) the bound is not violated. On the other hand if 
one considers the effects of higher orders, e.g. (9(A 8 T 4 ), since the corrections may decrease the 
lower bound of 0.257, there could be a possibility for violation of the bound in a narrow window 
where A 2 T <C 0.07. It would be interesting to explore this point better. 
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Appendix A 

The Fi(p; A) functions on the right hand side of equations of motion in (2.12) to (2.15) are given 
by 

p 1296 f~ 4 34 91273 r 6 M 2003147 ~ 4 101 ~ 4 30 9179675 ~ 8 19861 ~ 4 

Fl = —[ Xp + 15552 AP + ( ^3312" A " 108 )A P + ( ^79936" A + l37J A 

485 ^ 2 28 ,2742697-g 17461043 ~ 4 8569 ~ 4 26 3091 ~ 12 83656099~ 8 

~ 048 9 + ^ 139968 2239488 + 1944^ 9 + ~972~^ 7121664 

72409 ~ 4 4827,~ 2 24 5545 /Tlfi 116087~ 12 1530733~ 8 193517~ 4 
» A + Z7^) A V + T7^( A ~ inn „ n A + -————A - ———A 



148368 6182' r 648 v 19962 958176 133080 

432 22 206 ~ 16 132625~ 12 388271~ 8 2123~ 4 4075 ~ 2 20 14 63,ri 6 

+ 1109 )P + ^T (A + ^9328" A ~ 474624 A ~ 1648 A + 1648 )A P " 1944 (A 

1229863~ 12 5417857 TS 309279 T4 405, 18 407 /Tlfi 471325~ 12 

A + A 8 A 4 + )p 18 (A 16 A 12 

19152 93632 5852 77 IH 432 v 43956 

4332419~ 8 _ 1694084~ 4 13396 ~ 2 16 12197 ~ 12 44747~ 8 49288231 ~ 4 

" 1054944 32967 + 1221 ' 9 + 5832 ^ ~ 73182 + 1170912 

_ 1601631 - 4 14 137 - 12 122791 ~ 8 44678345 - 4 _ 411869 - 6 12 20891 - 8 

97576 ' 9 324 [ 13152 236736 3288 ' 9 8748 [ 

14499791~ 4 _ 141442003 ~ 8 10 15707 ~ 8 784154~ 4 _ 121617627 ~ 10 8 

+ 668512 5348096 ' P + 34992 ^ + 15707 1005248 ' 9 

409103 -4 _ 4562105 ~ 12 901 ( ~ 4 _ 3603757 ~ 14 _ 46409 - t _ 42715 - 18 

69984 1 818206 1 9 1296 1 194616 ' 9 15552 9 139968 



3 / ~ yo _ ( 933871 - 4 _ W_y x2p2 s _ ( 12450601 ~ 4 _ 365 )X y 6 + ( 57784183 ~ g 



5 V 621216 719' ' v 4969728 719' ' v 4969728 

38395 T4 24265 nT2 24 ,5721773 TR 287905 T4 230515, T4 22 108895, ~ 12 

+ A 4 )A 2 p 24 + ( A 8 A 4 )A 4 p 22 + (A 12 

6471 8628 ; H y 207072 138048 34512 ; H 5752 v 
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10478339-g 775829~ 4 61664 ~ 2 20 10465 ~ 16 582367~ 12 7134823~ 8 
" 23521320 ~ 980055 + 326685'' 9 + 1438 ^ ~ 347760 + 6027840 
4931 ~ 4 45 18 4425 r 16 38731 ~ 12 4828817 ~ 8 2516369 ~ 4 _ 647 - 2 16 
8372 2093 )P 2876 1 17700 424800 238950 4425 ' 9 

2790 ~ 12 409321-g 18377303~ 4 11389 ~ 4 u 341 ~ 12 9093~ 8 26462905~ 4 
+ 719 ^ + 66960 + 4821120 + 14880'' 9 + 719^ + 1352~ + 589248 

106413~ 6 12 ,62087 r 8 10224109 - 4 7861183 r 8 10 1301 ~ 8 562726 ~ 4 
+ T09lT )Ap + 17256 (A + 2235132 A + 1986784 )A 9 + 2876 (A + ^5127 A 

107584111 ~ 10 266509 ~ 4 11499791 ~ 126 4367 ~ 4 213407 - 14 4 
+ 2248128 )A 9 +( 155304 A + 621216 )A 9 + ( 25884 A + ^sW )A 9 
795431 r 16 , 93911 ~ 18 
310608 9 310608 



5 /~ 4 ^ 28 + 1751323 - 4 _ i0_ )pp26 + ( 22017197 ~ 4 + 403 + ( 988465 ~ 8 



3 V 481248 557' ' v 3849984 1671 y v 160416 

92237~ 4 1231,~ 222 80353 /T8 1926317~ 4 20513 N ~ 4 2n 8617 /Tl2 

+ A UV 2 + (A + A 4 )AV + (A 

80208 6684 ; F 13368 v 5785416 160706 ; P 2228 v 

2031685 T8 91487 ~ 4 2362 - 2 18 1587,r lfi 55925 Tl2 76121 ~ 8 

+ A 8 A 4 + )A 2 p + (A 16 + A 12 A 8 

1861272 155106 25851 7 F 2228 v 6348 16928 

25091~ 4 270, 1fi 5161. r 12 55141 ~ 8 663257 T4 7575 ,~ 2 u 463 /Tl2 

H A 4 )p 16 H (A 12 A 8 H A 4 UV 4 H (A 

9522 529 JH 167V 495456 743184 20644 ; H 557 v 

4315 T8 3655517 T4 25291 nT4 12 5821 ,~ H 125535 T4 223471 N r fi in 

+ A 8 + A 4 )A 4 p 12 + (A 8 + A 4 )A 6 p 10 

44448 800064 7408 ; ' 13368 v 23284 23284 ; H 

347 ,- s 99827 T4 16771699 N r 8 8 18701 ,~. 376574 N ~ in , , 301 T4 

+ (A 8 + A 4 )A 8 p 8 + (A 4 )A 10 p 6 + ( A 4 

2228 v 18738 599616 ; H 120312 v 18701 ' H v 20052 

_ 78325 ~ 12 4 _ 24203~ 14 2 _ 11795 ~ 16 

53472 ' 9 60156 9 240624 



= 324f AV 8 + — AV 6 + j 1893883^4 + Ifyy. + ( 3103^ 8 + 2747^ _ 20 - 
V 324 H y 559872 27 ; ' v 1458 1296 27 ; ' 

2479 ~ 8 121927 -4 _ 7135 ~ 4 20 1775 ~ 12 _ _6347_- 4 _ 125927~ 8 _ ^u2, 8 
648 1 535464 14874 ] 9 324 1 21300 383400 355 ' 9 
295 ~ 16 401 ~ 12 _ _6217_~8 _ 4199_r 4 _6_. 16 1139 ^ 2645 ~ 8 _ 14053 ~ 4 
108 1 1180 50976 5310 59 JP 324 1 4824 13668 
300 ~ 2 14 61 ~ 12 .lO^jg 520673 - 4 2819 ~ 4 12 1063 - 8 30295~ 4 

+ 1139 )AP + 81 (A + T756¥ A ~ 105408 A + W A P + "486~ (A " 51024 A 
H235 - 6 10 137 ~ 8 211 ~ 4 1529707 - 8 8 1921 - 129061 ~ 

+ ^504 )AP + 324 (A ~ 2466 A + "236736~ )A 9 + 17496 (A + ^684~ )A P 
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_79,r4 , 19495 ~ 12 4 3937 ~ 14 2 955 ~ 16 \ 
2916 1 632 ' P 17496 P 34992 J 



Appendix B 



A = 2(r 4 - r 4 )r + (^12073824A 2 r 18 + 21181440AV 8 ln(l + ^-) 

- 21181440A 2 r 2 r 16 - 21181440r 14 A 2 r 4 ln(2) - 2023875r 2 r 14 - 50400r 14 A 2 r V 
+ 27404704r 14 A 2 r 4 - 129 3 2 80r 12 A 2 r 6 + 39900r 10 A 2 r V + 307125r 6 r 10 

- 34650r 8 r 6 ^ 2 - 3686088r 8 r 10 A 2 + 14700r 6 A 2 r 10 ^ 2 + 5748750r 6 r 10 

- 19020720r 4 r 14 A 2 - 4032000r 2 r 14 + 17341800r 18 A 2 

B = ^(r 4 - r 4 )r + 1 ^°' 7 - ( 12073824A 2 r 18 + 21181440A 2 r 18 ln(l + 

2 lboOOr^ \ r z 

- 22970640A 2 r 2 r 16 - 21 181440r 14 A 2 r 4 ln(2) - 2023875r 2 r 14 + 201600r 14 A 2 r 2 ^ 2 
+ 27404704r 14 A 2 r 4 - 125 70 1 80r 12 A 2 r 6 + 336000r 10 A 2 r V + 838425r 6 r 10 

- 201600r 8 r 6 cj 2 - 5600588r 8 r 10 A 2 + 268800r 6 A 2 r 10 c; 2 + 5710950r 6 r 10 

- 19007420r 4 r 14 A 2 - 4088700r 2 r 14 + 17489500r 18 A 2 ^ 

C = 4 fsr 6 + A 2 r 4 r 4 - r 2 r 4 - A 2 r 8 ^) + ^ T . ( 2299500r 14 r 6 - 37910250r 10 r 
r z \ J 12600r i4 \ 

mmtinnnJ.. 18 , hrr,r.r»rr._8 10 2 , n^mrnnJ.. 14 riMfnnnrl. 18 \ 2 



10^, 10 

- 43218000r 2 r 18 + 790650r 8 r 1( V + 81742500r 6 r 14 - 546860235r 4 r 18 A 2 
+ 341269458r 8 r 14 A 2 - 493271 73r 12 r 10 A 2 - 50575808r 14 A 2 r 8 + 72793560r 16 A 2 r 6 
+ 20495744A 2 r V 8 - 122473890A 2 r 2 r 20 + 72442944A 2 r 22 + 262235400r 22 A 2 

2 

+ 127088640r 18 (r 4 - r 4 )A 2 ln(l + T -\) - 42362880r 14 r 4 (r 4 - r 4 )A 2 ln(2) 

- 22050r 12 r 6 cu 2 - 226800r 18 A 2 r 2 cu 2 - 2913750r 18 r 2 - 497700r 14 A 2 r 6 c 

- 1421700r 10 A 2 r 1( V - 1465800r 6 A 2 r 14 cj 2 ] , - 1 Ts 

J (r 4 - r 4 ) 
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D = - — : 1 J (-8 + AV)r 10 - (8A 2 r 4 + u 2 )r 8 + (8 - AV)r 4 r 6 + 8A 2 r 12 



2(r 4 - r 4 )r 



Ao^co^x2_i8/ 32782 4 2 6 100864 



100800^(r 4 - r 4 )W ^07382^ V 8 (-—, 4 .V - — r V + .V° 

+ iS r ° 8 ^ r2 + i^f r ° V + ill r ° 12) ln(1 + £ } + 42362880 ^o 8 r 20 m(2). 2 

+ 193181184A 2 r 2 (-|cj 2 + r 2 )r 26 + (-301985040A 2 r 6 + 72442944A 2 r V)r 24 

8 

+ (-216910848A 2 r 8 + 180332460A 2 r 6 a; 2 )r 22 + (588652400A 2 r 10 

- 151400000A 2 r 8 cj 2 + (4473000cj 2 - 756000A V)r 6 )r 20 + (23729664A 2 r 12 
+ (-4914000 + 9837212A 2 cj 2 )r 10 )r 18 + (1366241912A 2 r 14 + (-590625cj 2 

- 1314600AV)r 10 )r 16 + ((-45108000 + 9495348A V)r 14 + 466200r 1 V)r 14 
+ (-7974974584A 2 r 18 + (-331800AV - 15600375cu 2 )r 14 )r 12 + (48816276A 2 cj 2 
+ 159390000)r 18 r 10 + (6879600r 18 cj 2 + 12897439512A 2 r 22 )r 8 + (-163800000 

- 29432400A 2 cj 2 )r 22 r 6 - 8635726200A 2 r 26 r 4 + 54432000r 26 r 2 + 2060352000A 2 r 30 

E = ( 204A 2 r 16 - 317A 2 r 4 r 12 - lllr 4 r 10 + 299A 2 r 8 r 8 + 222r 8 r 6 

24 r ii y 

- 463A 2 r 12 r 4 - lllr 12 r 2 + 277A 2 r 16 

12r 12 

i2„ 124 iii 12 2 | n 7\2„ 16 



708A 2 r 16 - 753A 2 r 4 r 12 - 267r 4 r 10 + 805A 2 r 8 r 8 + 378r 8 r 6 



1047AV i V - llWV + 287AV J 



Appendix C 

Introducing the uo = 2tqV0 the coefficients in equation of motion are 

G x = -i7(p 4 -l)yfA 2 (p 8 + — p 4 + — )--p 2 
KH ' H V 204^ 204 ; 68^ 

C 2 = -^(p 4 - l) 3 p 2 (^1020A 2 p 12 + 703A 2 p 8 + lllp 6 - 194A 2 p 4 - 999p 2 + 3047A 2 ^) 
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= _25216 4 _ 2 / 1 isp ln(1+ 1 ) + I p i4p ln(2) _ ™! isp 
3 15 Kl J H \ 2 H v p 2 ' T y > 63040^ 

18053 T2 1fi ,19275 15 856397 2 ,, 14 100127 T2 12 

+ A 2 p 16 + ( 1 A 2 ( V to ))p H AV 2 

25216 r v 403456 197 v 100800 " F 302592 H 

505 ,3699 T2 2 , 1n , 15293 T2 555 2 , „ 445 ,3201 

+ ( h A 2 to 2 )p 10 + ( A 2 to 2 )p 8 H ( 

6304 v 6464 /H y 504320 12608 IH 6304 v 2848 

T2 2 , - 296643 4T2 35205 2 125365 T2 

+ AW)p 6 p 4 A 2 p 2 H A 2 

" 706048^ 100864^ 100864 

C = -^(p 4 -i) 2 (-ip 20 P(i + ^^ 

_ 17967 ~ 28141 ~ 3 1285 - 2367401 

63040^ 126080^ 197 l 2048 1 201600 ^ 

3 1285 ~ 2 2 10355621 14 303 4587 ~ 2 100127 

+ 197^2048 + ^° 806400 '' P " 6304^6464 + ^ 72720 

2 ,, 12 , 303 , 2 104993 N ~ 2 555 2 13761 , in 

+ to 2 ))p 12 -( (t» 2 )A 2 to 2 H )p 10 

^ v 6304 v 24240 ; 12608 403456 ;P 

623 2 14999, 555 9 10119, „ 623 , 10119 
-(A 2 (to 2 - — — ) - — -to 2 - —-^)p 8 



6304 v v 7120 ' 1246 2848 " 6304 v 2848 

T 2 , 5497533 2u fi , 32439 5497533 T2 , 4 , 149383 T9 

A 2 ( V t» ))p + ( A 2 )p 4 + ( A 

v 348880 " F y 100864 3530240 JH v 100864 

32439 , 2 149383 T2 



100864 /r 100864 

o/ 47912,2837 32782 4 s ,~ 2 2 i fil 1 , 47912 

C * = (5989 - 17967" < + , 8 )AW 6 ln(l + + — < 

6304 igr 2 2i /o\ i 22\2 2 ^ 20\2 2 143121 lg ~ 2 2 

-p A to m(2) + p A to — -p A to — p A to 



17967 r ' v ' r 5 r 95824 

i» (?| _5 + , v » v+l — ? ,. 

35934 v v 313 ; 10016 Jl V 3353840 

185 4 , 10 395 ,~ 2 , 62 4 , 46701 2 , 8 1675 ~ 2 
-to 4 )p 10 + (A 2 (-— + to 4 ) to 2 )p 8 A 2 



11978 /r 35934 v v 79 ' 2528 ,r 35934 

96259 6 ~ 2 2 3225 2 1480 ~ 2 4 86525 2 - 2 2 , 

+ p A to + ( to H A )p p A to 

3353840 P v 23956 17967 w 503076^ 
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